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Abstract. Numerical simulations are performed in the framework of nonlinear two 
- dimensional magnetohydrodynamics to investigate the influence of mass loading 
on the solar wind interaction with Venus. The principal physical features of the 
interaction of the solar wind with the atmosphere of Venus are presented. The 
formation of the bow shock, the magnetic barrier, and the magnetotail are some 
typical features of the interaction. The deceleration of the solar wind due to the 
mass loading near Venus is an additional feature. The effect of the mass loading is 
to push the shock farther outward from the planet. The influence of different values 
of the magnetic field strength on plasma evolution is considered. 


Introduction 

Data from the Pioneer Venus Orbiter (PVO), as well 
as data from Venera 9 and 10 and other sources, have 
contributed to a growing understanding of the inter- 
action of the solar wind with an unmagnetized planet 
such as Venus. The general picture derived from the 
available data is that, for typical solar wind conditions, 
an electrically conducting ionosphere deflects the super- 
sonic and superalfvenic solar wind around the planet. A 
bow shock forms upstream of the ionosphere and serves 
to slow and also assists in deflecting the solar wind. 
The size of the bow shock depends on the sonic and 
alfvenic Mach numbers and on the shape and the size 
of the effective obstacle, whereas asymmetries in the 
shock shape are determined by the direction of the in- 
terplanetary magnetic field (IMF) [ Khurana and Kivel- 
son , 1994]. Zhang et al. [1990] have proven that at 
Venus the effective obstacle undergoes a systematic size 
change of 0.13 R ^ (where R ^ is the radius of Venus) dur- 
ing a solar cycle. The boundary separating the shocked 
and magnetized solar wind plasma from the thermal 
ionosphere plasma is referred to as the ionopause, and 
the region between the bow shock and the ionosphere is 
referred to as the magnetosheath. The inner portion of 
the magnetosheath contains a region of enhanced mag- 
netic pressure referred to as the magnetic barrier [e.g., 

Luhmann, 1986]. There are physical processes involved 
in these relatively large-scale interaction regions as well 
as in the magnetic tail that the present work has appli- 
cation to. 
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Zhang et al. [1991] have used PVO data to show 
that the convected field gasdynamic model of Spreiter 
et al. [1966] predicts the correct bow shock location 
if the magnetic barrier is taken as the obstacle around 
which the plasma is deflected. They also found that the 
magnetic barrier is strongest and thinnest (about 200 
km thick) at the subsolar point and becomes weaker 
and thicker at the terminator. The magnetic barrier is 
responsible for transferring most of the solar wind dy- 
namic pressure to the ionosphere through the enhanced 
magnetic pressure. 

Observations at Venus have revealed that a magnetic 
tail forms as a consequence of the solar wind interac- 
tion with the Venusian ionosphere. The portions of IMF 
lines passing near the surface of Venus are slowed due 
to the interaction with newly ionized atmospheric neu- 
trals. The two ends of the magnetic ropes continue to be 
pulled downstream by the solar wind [e.g., Spreiter and 
Stahara, 1980]. Slavin et al. [1989] have estimated that 
the magnetotail extends to between 50 to 150 R^ from 
the center of the planet for a typical solar wind flow 
with the Alfven speed V& = 60 km/s and the Alfven 
Mach number Ma = 7.2. 

A study by Nagy et al . [1981] showed that Venus 
has a dayside neutral exosphere dominated by oxygen 
at altitudes above ~ 400 km calculated from the plane- 
tary surface. The exospheric oxygen together with other 
particles present can be ionized by solar radiation and 
charge exchange. Newly ionized particles are electro- 
magnetically coupled to the shocked solar plasma of the 
magnetosheath. As a consequence, near the ionopause 
the plasma becomes mass loaded due to interactions 
with ions, particularly O* [e.g., Luhmann et a/., 1991]. 
Model computations by Gombosi et al., [1980] using in- 
put from PVO ionosheath and neutral atmosphere pro- 
files indicate that approximately 10 % by number of the 
solar wind protons may undergo charge exchange with 
neutrals at ionopause altitudes. As a result of the mass 
loading (ML), the plasma is slowed and compressed. 
When the plasma then flows around the planet into the 
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magnetosheath, momentum conservation requires an in- 
crease in flow speed into the magnetotail. 

It has been noted by Luhmann [1986] that the above 
general behavior tends to prevail providing the thermal 
pressure in the upper ionosphere is sufficiently larger 
than the solar wind pressure and the pressure balance 
occurs where the ionospheric plasma is collisionless. As 
the solar wind pressure increases, relative to that of 
the ionospheric plasma, the magnetic barrier eventu- 
ally breaks down and the height of the dayside iono- 
sphere decreases. When the ionopause is too close 
to the planet, the transterminator transport of iono- 
spheric plasma ceases, thereby leading to a nightside 
phenomenon known as the disappearing of the iono- 
sphere [Cravens et a/., 1982]. 

Cloutier et al. [1987] have noted that ML processes 
at Venus should be asymmetric since the solar wind 
electric field depends on the relative directions of the 
flow velocity and the IMF. Phillips et al. [1987] have 
claimed that the flow and field configuration of the mag- 
netosheath plasma, together with the large gyroradius 
of the pickup ions, cause ML to occur preferentially on 
one side of the magnetosheath. More efficient pickup 
of newly created ions should occur over the hemisphere 
that produces an electric field directed outward, away 
from the ionosphere. Luhmann et al. [1991] and Phillips 
et al. [1986] have presented observations indicating an 
increase in the magnetic field intensity and current over 
the hemisphere where the electric field is directed out- 
ward. The enhanced solar UV during solar maximum 
increases the ionization rate of neutrals, which moves 
the bow shock out from the planet both in the equa- 
torial plane and at the terminator [Zhang et al ., 1990]. 
Similarly, Alexander and Russell [1985] have shown that 
the position of the bow shock at the terminator de- 
pends on solar activity. This dependence may be indi- 
rect evidence of ML since, when the solar EUV is high, 
more mass may be added to the shocked solar wind, 
thereby forcing the bow shock to recede from the planet. 
Linker et al. [1989] have revealed that ML can increase 
or decrease the plasma temperature, depending on the 
value of the sonic Mach number Ms . When 7 = 5/3, 
M s > a/ 9/5 is required for heating to occur. Numerical 
simulations performed show that the effects of ML on 
plasma temperature and velocity are most pronounced 
in a wake region. However, this theory was developed 
for the case of Io’s atmosphere, which is described by 
MHD equations with the source terms in the mass con- 
tinuity, Euler, and energy equations. In the case of Io, 
the magnetic field lines are in the north-south direction 
and are perpendicular to the horizontal flow. 

Recent computer simulations of the three-dimensional 
global interaction between the solar wind and unmag- 
netized bodies have proven to be extremely useful for 
improving our understanding of the associated large- 
scale processes [e.g., Wu, 1992; Tanaka , 1993; Cable 
and Steinolfson, 1995; Gombosi et al. 1994; McGary 
and Pontius , 1994]. Venus is of particular interest since 
the planet does not have a significant intrinsic mag- 
netic field, and the interaction of the solar wind with 
the Venusian ionosphere involves fundamentally differ- 


ent physical processes than occur at magnetized Earth. 
In addition, a large quantity of relevant observations are 
now available for Venus. When interpreted in associa- 
tion with the simulated results, this data set provides a 
test of the physical processes included in the model. 

Single-fluid MHD simulations of the global interac- 
tion of the solar wind with Venus without considera- 
tion of the ionosphere have been performed by several 
investigators. Wu [1992] limited his study to the day- 
side. Tanaka [1993] included the nightside as well but 
only considered magnetic field orientations parallel and 
perpendicular to the solar wind flow and was in a pa- 
rameter regime not directly applicable to average solar - 
wind conditions at Venus. Cable and Steinolfson [1995] 
carried out simulations for typical observed solar wind 
conditions at Venus and for an IMF orientation near 
the Parker spiral. McGary and Pontius [1994] studied 
the effects of ML, but they considered only the dayside. 
Moreover, they discussed the case of a cylinder with 
d/dz = 0 and B along the cylinder axis ( z - axis) and 
the perpendicular flow, namely, V 1 B. 

Our purpose is to extend the model of McGary and 
Pontius to study the nightside and to discuss various 
strengths of the magnetic field. We also considered a 
more realistic spherical model. 

The paper is organized as follows. The physical model 
used in the present study is described in the next sec- 
tion. Section 3 describes the numerical model which is 
used in the present studies. We present and discuss our 
detailed results in section 4. The paper concludes with 
a short summary. 

Physical Model 

We solve the following form of the compressible MHD 
equations as an initial value problem: 

+ V ■ (pV) = q 0 exp [-(h - h 0 )/ Ho), (1) 


d -^p- + v • [(pV)V] = - Vp + 2- (V X B) x B, (2) 


^ = Vx(VxB), 


( 3 ) 
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+ - [B x (V x B)]} = 0, 


( 4 ) 


V-B = 0. (5) 

Here p is the plasma density. V is the velocity. B is 
the magnetic field, p is the plasma pressure, and qo is 
the production rate at a reference altitude h 0 . H 0 is the 
scale height of the hot neutral oxygen population, and 
the ratio of specific heats is 7 = 5 / 3 . 
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In the present first-order approximation, the newly 
created photoions are considered as cold (T = 0) and 
motionless (V — 0) particles. A more rigorous approach 
would be to add corresponding source terms to the Eu- 
ler (2) and energy (4) equations [e.g., Gombosi et a/., 
1994]. However, it was shown by Biermann et al. [1967] 
that the major effect is that associated with the mass 
continuity equation and that contributions to the mo- 
mentum and energy density are small. Therefore in 
these preliminary calculations source terms in the mo- 
mentum and energy equations are neglected. The ML 
term is given by the right-hand side of the mass continu- 
ity equation (1). Oxygen photoionization is considered 
to be the only mass source with an altitude distribution 
described by the exponential function. This production 
rate is based on PVO observations for solar maximum, 
where the hot neutral oxygen atmosphere had a scale 
height Ho = 400 km in the subsolar region and a refer- 
ence production rate of go — 3 x 10 5 cm -3 s“ 1 at altitude 
ho = 400 km [Belotserkovskii et a/., 1987]. As numerical 
simulations, which were made for various mass addition 
rates, revealed that the boundary layer develops when 
oxygen ion production rate is q 0 = 3 x 10 5 cm~ 3 s _1 
[McGary, 1993; McGary and Pontius, 1994], we carried 
on our simulations for this value of q 0 . 

Numerical Model 

These numerical calculations were performed with 
HEMIS3D, a three - dimensional ideal MHD code. The 
code utilizes the modified Lax-Wendroff scheme devel- 
oped by Rubin and Burstein [1967], which yields stable 
solutions by adding numerical diffusion to the scheme. 
The dissipative terms (kept as small as possible) are 
included to help damp out short-wavelength ripples 
generated by numerical dispersion and numerically in- 
duced reflections from the simulation boundaries, while 
leaving the longer- wavelength phenomena minimally af- 
fected. The equations are solved in all three spatial di- 
mensions of a spherical ( r , 0 , <f>) coordinate system in 
which the 0 = 0 axis is directed toward the Sun or into 
the solar wind flow. The generally small deviations of 
the solar wind flow from the radial direction (with re- 
spect to the Sun) are neglected so the solar wind flow 
impacting the Venusian ionosphere can be assumed par- 
allel to the Sun- Venus line (8 — 0). The spherical co- 
ordinate system is centered on the planet so r = R v 
represents the planetary surface. 

Although a grid defined in spherical coordinates in- 
troduces some computational complexity, it has a num- 
ber of advantages for studying flow past planetary ob- 
stacles such as Venus. One obvious advantage is that 
the boundary at the planetary surface occurs at a fixed 
value of the radial coordinate, which avoids the compli- 
cations of interpolating in a Cartesian coordinate sys- 
tem. Another distinct advantage is that as the center 
of Venus is located at the center of the spherical coor- 
dinate system, high resolution near Venus is obtained 
without the expense of high resolution everywhere in 
the simulation. The code is constructed in such a man- 


ner that it can be applied to a single r, 0 plane when 
the phenomenon being studied is cylindrically symmet- 
ric about the poles at 0 = 0 and 8 = 180°, as is the case 
when the magnetic field is parallel to the solar wind 





Figure 1. Plots of the normalized mass density profiles along 0 
= const lines for p = 0.1, p = 0.5, p =10 and for the NL and 
ML simulations, (a) the Sun-Venus (0 = 0) line, (b) the 
terminator (0 = 90°) line, and (c) the 0 = 180° line. Mass 
loading shifts the bow shock farther away from Venus. 








MURAWSKI AND STEINOLFSON: SIMULATIONS OF MASS LOADING 


2551 


Plate 1. Normalized mass densities (p/poo) for the MHD non loaded (top) and mass loaded 
(bottom) simulations for the case of a parallel (B||V) magnetic field and the plasma (3 = 0.5. 
Here, is the density of the solar wind. The effect of mass loading (ML) is to increase the 
mass density in the overall region, except the tail region in which the mass density is reduced. 
Poo denotes the solar wind mass density. 


CFL condition, At < A//|V|, where the cell size is A / 
and V is the fastest physical speed in the system. The 
resulting computer code [Cable and Steinolfson, 1995] 
is second-order accurate in space and time. 

Defining a spherical boundary for the problem is 
easily accomplished. However, solving the equations 
of magnetohydrodynamics in spherical coordinates re- 
quires some modification of the usual Lax-Wendroff 
scheme. In particular, the singular points at 0 = 
0 and 0 = 180° are handled through application of 
L’Hospital’s rule to extrapolate values from grid points 
next to the singular points. Since the value of each vari- 
able at the singular points may vary depending on the 
angle <f> at which the singular points are approached, a 
flux-weighted average of the values at all angles is used 
for the final updated value. This procedure is described 


in more detail by Cable and Steinolfson [1995]. 

We used a simple model for the boundary conditions 
that is physically motivated and performs adequately 
well. Venus is modeled as a hard, highly electrically 
conducting sphere. The inner boundary is taken at the 
planetary surface (1 i^) where the radial components of 
both the magnetic field and velocity are set to zero. All 
other physical variables are computed by extrapolation 
of values along a radial line. At the day side of the outer 
radial boundary, all quantities are set to the solar wind 
values. The values at the nightside of the outer radial 
boundary are determined by extrapolation along the 
local flow direction. Moreover, we have never encoun- 
tered transient conditions where the flow was coming 
in from the outer radial boundary. A typical compu- 
tation begins with the introduction of the desired solar 




2552 


MURAWSKI AND STEINOLFSON: SIMULATIONS OF MASS LOADING 


Plate 2. Normalized plasma flows ( V/a <») for the MHD non loaded (top) and mass loaded 
(bottom) simulations for the case of a parallel (B||V) magnetic field and the plasma f3 — 0.5. 

Here, is the sound speed of the solar wind. The effect of ML is to decrease the plasma flow. 

wind values in the dayside within the numerical box. is initialized in a similar way, except that B is skewed 

The initial IMF is taken from the potential solution for with respect to the solar wind flow by the desired angle, 

flow over a sphere so that V ■ B = 0 in the initial state. In the present paper, this angle is 0°. This is clearly 

The continued satisfaction of V • B = 0 is assured by not a frequently occurring case [Phillips et a/., 1986], 

solving a Poisson equation and updating the magnetic but it is a reasonable first-order approximation which 
field. The numerical solution then continues until the will serve as a reference for comparison with results for 

interaction achieves an approximate steady state. cone angles of 45* and 90°. These processes will be 

The density and pressure are initially constant through- realistically modeled in future calculations, 
out the simulation region. The initial velocity field is 
the flow of an incompressible fluid over a sphere, di- 

rected along the d = 0 line: Results and Discussion 

V r (r,6,<f>,t = 0) = -Vb[l — ( — ) 3 ]cos0, (6) We Present all numerical results for the oxygen ion 

r production rate q 0 — 3 x 10 5 cm~ 3 s _1 [McGary and 

1 ry Pontius, 1994] and discuss several values of the plasma 

Vg(r, 6,<l>,t= 0) = Vb[l + £(— ) 3 ]sin0, (7) 0 = Sirpo/Bl namely 0 = 0.1, 0 = 0.5, 0 = 10, and a 

r typical value of 0 = 1.5. The case of 0 = 0.1 (/? = 10) 

corresponds to the strong (weak) magnetic field case. 
The physical solar wind values used are: electron den- 
where Vq is the solar wind speed. The magnetic field sity n t = 20 cm -3 , temperature T = 10 5 °K, sound 


= o) = o, (8) 
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Plate 3. Normalized plasma pressure {p/{p 00 a 2 00 )) for the MHD non loaded (top) and mass 
loaded (bottom) simulations for the case of a parallel (B||V) magnetic field and the plasma 
fl = 0.5. The effect of ML is to decrease the plasma pressure. 


speed 55 km/s, flow velocity 370 km/s, sonic Mach num- 
ber 6.73, and thermal pressure 6.06 • 10" 10 dyn/cm 2 . 
Other physical parameters are chosen the same as in 
the work by Cable and Steinolfson [1995]. 

As a first step to apply the methodology to the prob- 
lem of ML in the atmosphere of Venus, we consider an 
axisymmetric case with the magnetic field lines parallel 
to the solar wind flow. Consequently, B is parallel to 
the Sun- Venus line and the entire solution is rotation- 
ally symmetric around this line, 8/d<j) — 0. We used 
the numerical technique described in the previous sec- 
tion to obtain a solution of the set of (l)-(5). The flow 
is from left to right and all physical quantities are either 
symmetric or antisymmetric about the lower boundary. 

As shown in Plate 1, the mass density is increased in 
the case of mass loaded (ML) solutions. This finding 
is consistent with a simple analysis of mass continu- 
ity equation (1) and with previously published results 


[e.g., Linker et al, 1989; McGary and Pontius, 1994]. 
However, a more detailed close-up reveals that in the 
tail and terminator regions the mass density is reduced 
by the ML effects (Figure 1). Similar reduction is ob- 
served in the case of (3 = 0.5 along the Sun- Venus line, 
at r ~ 1.07 R v (Figure la). The density reduction 
does not occur at the terminator for (3 = 0.1, while 
this effect takes place for ft = 0.5 and / 3 — 10 (Fig- 
ure lb). The density depletion is stronger for a higher 
value of (3. Just at the planetary surface, at r < 1.3 
itp, the mass density increase is observed (Figure lb). 
In the magnetotail, the plasma density reduction occurs 
at the planetary surface for all values of the plasma /3. 
Again, the rate of density depletion is proportional to 
(3. The case of (3 = 0.1 differs from the other cases just 
as at the planetary surface p is increased, but farther 
on there is a region of mass depletion which is followed 
by a region of mass enhancement (Figure lc). To our 
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Plate 4, Normalized plasma temperature (T/T 0 o) for the MHD non loaded (top) and mass 
loaded (bottom) simulations for the case of a parallel (B||V) magnetic field and the plasma 
(3 = 0.5. The effect of ML is to decrease the plasma temperature. 


knowledge, this is the first report of mass depletion due 
to the ML. The density depletion is a consequence of 
the nonlinear term V • (pV), which is smaller due to 
the plasma flow reduction by the ML effects. Indeed, 
the largest density (velocity) increase (decrease) occurs 
in the magentotail at r = 3 R v for (3 = 0,1. The nor- 
malized density jumps from pjvl — 1-25 to pml — 1.5 
(Figure lc), while the magnitude of plasma flow is re- 
duced from V^l — 4.9 to Vml — 4 (Figure 3c). Con- 
sequently, pnlVnl — 6.125 > pmlVml — 6. Therefore 
the effect associated with this nonlinear term overcomes 
the density increasing effect exerted by the exponential 
term goexp [— (k — h 0 )/Ho\. The density reduction in 
the tail may also be due to the changed shape of the 
bow shock. A different bow shock shape may result in 
less mass being carried into the tail, regardless of other 
factors. 

It follows from Euler equation (2) that we can expect 


the reduction of a velocity magnitude because conserva- 
tion of momentum and energy requires that the plasma 
flow is decelerated as a consequence of the mass density 
reduction. Plate 2 presents results both for the ML and 
no loading (NL) simulations. By examining the two so- 
lutions, one can see a fundamental difference between 
them. The effect of ML is to decrease the flow. A de- 
creasing rate is highest in the tail region where slow 
flows in the case of NL are much more reduced by ML. 
The decreasing rate is larger for a smaller value of (3 at 
the terminator and in the tail region. The plasma starts 
decelerating because of ML well ahead of the shock. In 
the subsolar region, this rate is largest for (3 — 0.5. A 
large flow deceleration occurs also at the terminator, 
close to the planetary surface where the plasma flow 
drops from V NL ~ 5.7 to V ML ~ 4.2 (Figure 2b). In- 
spection of Plate 2 reveals the global structure of the 
shock and the stagnation region behind it. The flow 
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decelerates at the shock as approaching solar wind ab- 
sorbs an increasing number of ions. At the planetary 
surface, r = iJ,, there is no perpendicular flow as a con- 
sequence of the boundary conditions imposed. Velocity 
shears in the radial direction are present at the bow 
shock region both for the ML and NL flows, which is 


consistent with PVO observations [Mihalov et ai } 1982] 
and the gasdynamic results by McGary [1993]. 

From (4) we can derive an equation which governs 
the evolution of the NL plasma pressure p 

^ + V • (pV) = (1 - y)pV ■ V. (9) 



Figure 4. As in Figure 1 but for the normalized temperature Figure 5. As in Figure 1 but for the normalized magnitude of the 
profiles. magnetic field profiles. 
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As the plasma flow velocity V is decreased by the ML 
effects we can expect a pressure reduction. The simula- 
tion results shown in Plate 3 are qualitatively consistent 
with our expectations. However, Plate 3 shows that the 
plasma pressure is increased near the terminator region, 
in the region 0 < 90°. This increase is caused by a 
nonlinear interaction which is described by the V * (pV) 
and pV • V terms. Our results differ from the conclusion 
made by McGary and Pontius [1994], who claimed that 
the dynamics of the plasma flow are not significantly 
affected by ML. A reason for the difference between the 
present results and those of McGary and Pontius has 
to do with the different models applied; McGary and 
Pontius performed gasdynamic (with B = 0) simula- 
tions for the cylinder flow, whereas our model concerns 
MHD and a plasma sphere. 

The above scenario may differ for different values of 
(3. It occurs in the case of a sufficiently small value of (3 
that the magnetic field is so strong that the fast shock 
becomes an intermediate one [Steinolfson and Hund- 
hausen, 1990]. For example, the ML increases the gas 
pressure along the whole Sun-Venus line in the case of 
(3 = 0.1 and (3 = 10 (Figure 3a) but reduces it in the 
case of (3 = 0.5, f3 = 1, and (3 = 1.5 at the planetary 
surface. In the bow shock region, an enhancement of 
p due to the ML is observed for all values of (3. The 
strongest enhancement seems to occur for (3 — 0.5. At 
the terminator (Figure 3b), the gas pressure is increased 
by ML. The rate of decrease is highest at the planetary 
surface for (3 — 0.1, while at larger distances the influ- 
ence of ML on p is negligible. The pressure reduction at 
the planetary surface is caused by the right-hand side 
term of (9) which becomes small at the surface because 
of the flow reduction. The bow shock is closer to the 
planetary surface for a higher value of /?, in accordance 
with our expectations. At the magnetotail the pressure 
is reduced (Figure 3c). The reduction rate is highest at 
the planetary surface for the smallest value of (3 and is 
a monotonic function of j3. 

The next point of comparison between the ML and 
NL results is the temperature distribution around Venus. 
As for the ideal gas the plasma temperature T is pro- 
portional to p/p so both p and p play a major role in 
determining its behavior. We already learned that for 
(3 — 0.5 p (p) is increased (decreased) by the ML effects 
everywhere except in the terminator and tail (nose) re- 
gions. Therefore, we may expect the general behavior of 
T would be to decrease its magnitude by the ML effects. 
Indeed, Plate 4 indicates that the ML significantly cools 
the flow relative to the NL case. The cooling occurs as 
zero-temperature oxygen ions are introduced into the 
flow and thermalized [McGary and Pontius , 1994]. The 
behavior of T is interesting to discuss in more detail as 
a function of the radial direction r at different locations 
of 0. In the case of 0 — 0, for (3 = 0.1 and f3 = 10 
the temperature is decreased by ML both at the plane- 
tary surface and at the sunward side of the bow shock 
(Figure 4a), while in the shock the effect of ML is neg- 
ligible. The curves which correspond to the ML possess 
maxima in the magnetosheath, while no maxima exist 
for the NL case. A similar maximum occurs also in the 
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Figure 6, Terminator positions for different values of the plasma 
p and for the non loaded and mass loaded simulations. The bow 
shock at the terminator is shifter farther away farther away from 
the surface of Venus. 

case of (3 — 0.5 for the ML curve. Now, however, ML 
cools the plasma at the planetary surface but warms it 
in the bow shock region. The temperature increase at 
the bow shock was observed recently by Gombosi et al. 
[1994] in their simulations of the interaction of an ex- 
panding cometary atmosphere with the solar wind. At 
the terminator, the plasma is warmed by ML (Figure 
4b). The warming is more significant at the planetary 
surface and is higher for a higher value of /3. The cooling 
occurs at the magnetotail (Figure 4c). The cooling rate 
is higher at the surface of Venus and for a smaller value 
of (3. This is a consequence of the pressure decrease 
at the planetary surface (Figure 3c). The mass den- 
sity is increased for (3 = 0.1 and decreased for (3 = 0.5 
and (3 — 10, but the decreasing rate is smaller than 
the pressure decreasing rate (Figure lc). Consequently, 

T ^ p/p is much decreased at the surface of Venus. The 
cooling rate is much more significant than the warming 
effect. The bow shock position is closest to the plane- 
tary surface for the smallest presented values of (3. The 
bow shock distance from Venus grows with (3 . 

Plate 5 shows the magnitude of the magnetic field. 
In the subsolar region there is no jump in the magni- 
tude of the magnetic field across the bow shock. At the 
flanks, however, the magnetic field increases through 
the bow shock. The magnetic tail is a region of reduced 
magnetic field. It can be seen that the ML exerts some 
influence on the magnetic field strength. Small differ- 
ences occur at the sunward side of the shock where the 
magnetic field is decreased by the ML effects (Figure 
5a). This effect is small for (3 — 0.1 and grows with 
(3. The decrease in the magnetic field strength nearest 
the planet along the stagnation streamline, is the re- 
sult of the boundary conditions (i? r (r = i2y,0, t) = 0) 
imposed together with the fact that the parallel shock 
does not change a value of the magnetic field across 
the shock. As the planet is conducting there can be 
no radial magnetic field component at the surface of 
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Venus. But near the subsolar line, the solar wind mag- 
netic field is purely radial. Consequently, the magnetic 
field strength has to drop to zero in this region. This 
behavior counter to the common picture of B piling up 
in the stagnation region in the case of B 1 V [Cable 
and Steinolfson, 1995; Murawski and Steinolfson, 1995] 
and would only occur for the magnetic field nearly par- 
allel to V. At the terminator, the ML enhances the 
magnitude of B at the bow shock. Therefore at the 
planetary surface B is reduced. A decreasing rate is 
hardly dependent on (3. The largest differences occur in 
the magnetotail where the magnetic field is enhanced 
by ML. The strongest enhancement occurs for (3 — 0.1. 
This behavior contradicts our expectations, which fol- 
low a simple analysis of (3) wherein as a consequence of 
V decrease we expect the magnetic field reduction by 
the ML effects. However, V couples nonlinearly with 
the magnetic field. The coupling occurs through the 
term V x (V x B), which determines the overall behav- 
ior of the magnetic field rather than the plasma flow V 
only. 

Inspection of Plates 1-5 reveals that a bow shock is 
formed upstream of Venus. The principal effects are 
that the bow shock is displaced farther from Venus in 
the case of the ML plasma [Belotserkovskii et a/., 1987; 
Breus et al. t 1989; McGary and Pontius , 1994]. For 
precise investigation of the bow shock structure, Fig- 
ures 1-5 exhibit the plasma mass density, velocity, gas 
pressure, temperature, and magnetic field distribution 
along the Sun- Venus (0 — 0°), terminator (0 — 90°), 
and 0 = 180° line. Inspection of these figures reveals 
that the pressure at the planetary surface is larger for a 
smaller value of/?. The ML effects also shifts the shock 
at the terminator farther away from the planetary sur- 
face (Figure 6). This shift is largest for intermediate 
values of the plasma /?. These results show the bow 
shock at the terminator is in about the position of the 
PVO observations [Tatrallyay et a/., 1984; Russell et 
a/., 1990]. Average subsolar and terminator distances 
of the bow shock are about 1.15 R v and 2 R v (from 
the center of the planet), respectively. To match the 
observed locations (at the nose) of the bow shock, the 
MHD model requires that the mass loading be greatly 
enhanced. 

Russell and Zhang [1992] in their analysis of PVO 
data found that during periods when the solar wind 
magnetosonic Mach number is near unity and the plasma 
beta is low, the bow shock may travel to distances 
greater than 10 R v away from the planet. On the ba- 
sis of previous work on the formation of MHD shocks 
in coronal mass ejections [Steinolfson and Hundhausen, 
1990], it was felt that the solar wind conditions had be- 
come (during the far bow shock excursions) such that 
the usual fast MHD shock is no longer a stable solu- 
tion. The lower than normal plasma beta and mag- 
netosonic Mach number are in a parameter regime for 
which the usual fast-mode bow shock close to the planet 
may not provide the necessary compression and deflec- 
tion of the solar wind. Using MHD simulations, it was 
shown [Steinolfson and Cable , 1993] that, for these con- 
ditions, the usual fast shock is replaced by a bow shock 
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Figure 7. Time history of the radial velocity component (V r /a°°) 
at the spatial position very close to the bow shock for the non 
loaded simulations. The upper (lower) curve corresponds to P = 
10 (p =0.1). Not the stabilizing effect of the stronger magnetic 
field. 

configuration containing an intermediate shock near the 
Sun- Venus line and a fast shock at large distances from 
the Sun- Venus line. The resulting configuration propa- 
gates upstream away from the planet at a low speed and 
appears to be approaching a new equilibrium stand-off 
location at a large distance from the planet. 

Plates 1-5 exhibit waves which propagate along the 
bow shock boundaries. We have checked that these 
waves are not numerical artifacts as they do not dis- 
appear for a finer numerical grid and are persistent in 
time (not shown). These waves are generated by the 
solar wind impact, which is exerted by the solar wind 
on the plasma surrounding the planet. The first impact 
shifts the plasma towards the surface of Venus, while 
later in time the magnetic field recovers and the plasma 
is displaced farther away from the planet. Moreover, 
the bow shock first builds closer to the planet and then 
its altitudes at the equator (0 — 0°) and at the termi- 
nator (0 = 90°) grow in time. As a consequence, the 
nose of the bow shock oscillates in time. The Alfven 
and slow oscillations, which are guided by the magnetic 
field lines [e.g., Murawski et a/., 1995], are fed by the 
Kelvin-Helmholtz (KH) instabilities [e.g., Rankin et al. } 
1993], which are more apparent in the case of a weak 
magnetic field case. Figure 7 presents the radial com- 
ponent of the plasma flow as a function of time and for 
(3 = 0.1 and /? = 10. This component is calculated at 
a spatial position which is very close to the bow shock. 
As expected, the stronger magnetic field makes the flow 
less unstable as the amplitude of oscillations is smaller 
in the case of (3 — 0.1 than in the case of (3 — 10. The 
KH instabilities criterion [e.g., Rankin et a/., 1993] 
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Plate 5. Normalized magnitude of the magnetic field (B / yj popoo^lo) for the MHD non loaded 
(top) and mass loaded (bottom) simulations for the case of a parallel (B||V) magnetic field and 
the plasma (3 — 0.5. The effect of ML is to increase the magnitude of B in the terminator region. 


is satisfied easier for a less magnetized flow. In this for- 
mula indices 1 and 2 correspond to variables on different 
sides of the bow shock, while denotes the growth rate 
of the KH instabilities. It is easy to find out that the 
flow is unstable (wf > 0) for a non- magnetic plasma 
(B X - B 2 = 0). On the other hand, strong magnetic 
fields can make the flow stable (a;? < 0). Intermediate 
values of the magnetic fields stabilize the flow and in 
particular the flow can be marginally stable (w? — 0). 
These instabilities occur in the region where the flow 
has a large gradient such as at the bow shock. Assum- 
ing that the wave propagation k is parallel to the flow 
V and magnetic field B, it follows from (10) that 


the units of the solar wind. Va is the Alfven speed of 
the solar wind. At the surface of Venus, V 2 = #2 = 
0. Upstream of the bow shock, the mass density pi 
and the plasma flow Vi are hardly dependent on the 
plasma f3. Therefore we can put p\ ~ 0.95 and Vi ~ 6.7 
in normalized units. For (3 — 0.1 (/ 3 = 10) we have 
p 2 — 3.9 (p 2 ~ 3.7) and B\ ^ 0.93 (B\ ~ 0.91). It is 
easy to check that the right-hand side of (11) is positive, 
indicating that the flow is KH unstable for both (3 = 0.1 
and j3 = 10, 

Summary 


The interaction of the solar wind with Venus is stud- 
2 o>? 2 2 2 2 ied using numerical solutions of the time-dependent 

(P 1 +P 2 ) — Pi P2 ( U 2 — Vx ) — (pi +p 2 )\B l -f B 2 ), ideal two - dimensional MHD equations. A modified 

(11) Lax-Wendroff scheme [Cable and Steinolfson, 1995] is 
where now the plasma parameters are normalized by used to solve the equations in a spherical coordinate 
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system. For these computations the detailed chemistry 
of the ionosphere is neglected, and the planet is treated 
as a conducting sphere. Numerical computations have 
been perfomed to study the effect of ML and the IMF 
strength on the magnetic barrier and the general con- 
figuration of the magnetic tail. The IMF orientation 
included in this study is parallel to the solar wind flow 
direction. 

The main results are the following: The ML increases 
the mass density in the overall region, except the termi- 
nator and tail where the mass density depletion occurs 
as a consequence of the nonlinearity action. The mass 
density depletion is larger for a higher j3. As a conse- 
qence of this reduction, the plasma flow, pressure, and 
temperature are reduced. This scenario depends, how- 
ever, on a value of the plasma (3, It was found that for 
/? = 0.1 both the plasma pressure and the temperature 
are increased at the day side by the ML. 
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